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ABSTRACT 

Heat-transfer through weakly magnetized diffuse astrophysical plasmas excites 
whistlers. This leads to electron whistler resonant scattering, a reduction of the elec- 
tron mean-free path, and heat-flux inhibition. However, only whistlers propagating at 
a finite angle to the magnetic field (off-axis) can scatter the heat-flux carrying elec- 
trons. Thus, the level of heat flux-inhibition along the magnetic field lines depends on 
the presence of off-axis whistlers. 

We obtain a solution of the Boltzmann equation with the whistler wave equation 
and show that if e th (3 e 3> 10~ 4 , where e th is the thermal Knudsen number, and j3 e is 
the ratio of the electron pressure to the magnetic energy density, scattering of heat- 
flux carrying electrons by off-axis whistlers, which are shown to propagate at about 
65°, is efficient enough to lead to heat-flux inhibition along field lines. The inhibition 
so obtained is proportional to (e (3 e )~ l . 

Key words: Heat-conduction, Whistlers 
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1 INTRODUCTION 

The presence of a cold gas embedded in a hot gas is a com- 
mon phenomenon. This situation occurs in many astrophys- 
ical systems such as the interstellar medium (ISM), the in- 
tracluster gas in galaxy clusters (ICM), AGNs and Lyman 
a systems. These multi temperature-density configurations 
were and still are the subject of intensive theoretical study. 
Field (1965) carried out a linear stability analysis of opti- 
cally thin, collision dominated magnetized plasma. He found 
the critical wavelength (the Field length) that distinguishes 
perturbation wavelengths that form cold clumps from those 
that don't. Perturbation wavelengths shorter than the Field 
length are heated by heat-conduction and fail to form cold 
clumps whereas the longer perturbation wavelengths are 
cooled enough by radiation losses to condense into cold 
clumps. The Field length retains its qualitative role in the 
non-linear regime (Cowie & Songaila 1977; Cowie & McKee 
1977; McKee & Cowie 1977). 

At large Knudsen numbers there is heat-flux satura- 
tion. This is considered by Cowie & McKee (1977), Bal- 
bus & McKee (1982), Slavin & Cox (1992), and Dalton & 
Balbus (1993). Chun & Rosner (1993) consider the effects 
of a non-local Maxwellian electron distribution on Field's 



analysis, and Bandiera & Chan (1994a,b) generalize their 
work and emphasize the role of the non-local thermoelec- 
tric effect in thermal evaporation processes. Apart from the 
details considered in these works, the Field length retains 
its qualitative and quantitative role. All the collision dom- 
inated plasma theories predict that scales smaller than the 
Field length evaporate (heated). For typical ISM or ICM 
parameters, the thermal evaporation theory predicts a hot- 
ter ISM and ICM than what is actually observed. This is 
also the case for Lyman-a systems and the standard AGN 
model. 

For the ISM, McKee & Ostriker (1977) proposed that 
the observed cold gas is constantly replenished by super- 
novae explosions. However, this does not necessarily apply 
to other astrophysical environments in which cold and hot 
gas are observed to co-exist. In broad line emission regions 
surrounding quasistellar objects (Begeleman & McKee 1990; 
McKee & Begeleman 1990) cold and hot gas coexist, but the 
presence of a small scale replenishment source is not clear. 
The ICM at the cores of some clusters is composed of inter- 
acting cold and hot gas (Canizares, Market & Donahue 1988; 
Fukazawa et al. 1994). Detailed studies show that thermal 
evaporation as predicted by classical theory (Spitzer 1962) 
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should be very efficient in the ICM ( Cowie & Binney 1977; 
Binney & Cowie 1981; Fabian & Nulsen 1977; Loewenstein & 
Fabian 1990; Fabian 1994). These works interpret the ICM 
observations as a reduction of the magnitude of the thermal 
conductivity below that predicted by classical theory. The 
evidence for heat flux inhibition at the ISM, ICM, AGNs, 
and Lyman-a systems is indirect. Direct conformation of 
heat flux inhibition exists in the solar corona and wind for 
three decades (Montgomery, Bame & Hundhausen 1968). 
Recent measurements performed with Ulysses (Scime et al. 
1994) confirm that between 1.2 and 5.4 AU, heat flux is 
inhibited compared to its Spitzer value. 

This work concentrates on the ICM observations, with 
a strong emphasis on the "cooling flow model" (CFM; cf. 
Fabian et al. 1984, 1991 Fabian 1994 for a review). The CFM 
reproduces the observed properties of the ICM in the central 
part of galaxy clusters very well, but only with the assump- 
tion of strong heat flux inhibition, i.e. heat flux is ignored. 
The justification of this hypothesis has been heuristic, and 
has presumed that heat flux inhibition is obtained by a tan- 
gled magnetic field ("in a manner which is not yet clear"). 
This paper considers heat flux inhibition in the frame work 
of a rigorous discussion, and attempts to quantify the level 
of heat flux inhibition. In doing so we partially rely on works 
that consider heat flux inhibition in the solar wind ( Gary 
& Feldman 1977; Gary et al. 1994 and Scime et al. 1994), 
due to whistlers. 

The ICM observations require heat conductivity sup- 
pression (heat flux inhibition) compared to the Spitzer value 
by at least two orders of magnitude, and sometimes up to 
four orders of magnitude (Pistinner & Shaviv 1996; Pistin- 
ner, Levinson & Eichler 1996). 

Reducing the electron mean free path would lead to heat 
flux inhibition. Mechanisms that would accomplish this can 
be divided into two main theoretical classes: 

(i) Increase of the path length of the magnetic field lines, 
hence the length over which parallel heat transport must 
operate. 

(ii) Plasma waves that resonantly scatter the heat con- 
ducting electrons. 

The first mechanism is realized when a tangled magnetic 
field pervades the plasma (for details cf. Cowie & Binney 
1981, Rosner & Tucker 1989; Tao 1995; Pistinner & Sha- 
viv 1996). The second comes about when self excited elec- 
tron plasma "waves exist" (Gary & Feldman 1977; Jafelice 
1992; Levinson & Eichler 1992; Pistinner et al. 1996). These 
plasma waves are expected to scatter the heat carrying elec- 
trons, thereby reducing their mean free path. 

The classification made above is quantitative (Gold- 
stein, Klimas & Sandri 1975), and it follows from the dif- 
ferent physical processes involved. More specifically, let Ib 
and ri be the coherence length of the magnetic field and the 
electron gyroradius respectively, then: 

(i) Mechanism (i) (the tangling mechanism) is realized 
when the condition rj <C Ib holds, and does not scatter 
electrons (except for adiabatic mirroring). 

(ii) Mechanism (ii), in which electromagnetic plasma 
waves scatter electrons "non-adiabatically" works if rt w Ib- 

Large scale tangling, (i), is now seriously constrained 
and in any case "questionable" (Tao 1995; Pistinner & Sha- 



viv 1996). It seems that only a very particular family of 
magnetic field correlation functions can inhibit heat con- 
duction. These correlation functions have Lorentzian rather 
than Gaussian properties, and the physical basis for their 
origin is not yet clear. In contrast electromagnetic elec- 
tron plasma waves (ii), obtained from a collisionless plasma 
model (Levinson & Eichler 1992), provide large inhibition 
factors when applied to the CFM (Pistinner, et al. 1996). 
The model of Levinson & Eichler (1992) follows Gary & 
Feldman (1977) and invokes whistler-electron interaction to 
suppress the heat transfer. In light of our findings (Pistin- 
ner et al. 1996) and the observational information and inter- 
pretation accumulated during the past three decades, heat 
conduction theory in diffuse astrophysical plasmas should 
be reexamined. Scattering by plasma turbulence should be 
considered along with Coulomb collisions. 

Plasma turbulence results from collective phenomena, 
and may lead to electron-waves resonant pitch-angle scat- 
tering. Non-magnetized plasmas are studied in detail, and 
under some conditions ion acoustic turbulence leads to heat- 
flux inhibition (Galeev & Natanzon 1984). Jaffliche (1992) 
considers this possibility in the context of the CFM. Al- 
though he ignores strict requirements for ion-acoustic waves 
excitation, he finds insufficient heat inhibition. Levinson & 
Eichler (1992) follow Gary & Feldman (1977) and consider 
heat flux inhibition by whistlers. Gary & Feldman (1977) 
find that the inhibition obtained by whistler-electron scat- 
tering in the (significantly magnetized) solar wind is not 
large. Levinson & Eichler (1992) show that collision domi- 
nated plasmas are subject to Weibel type instabilities that 
generate the whistler electromagnetic mode effectively, and 
find strong heat flux inhibition in weakly magnetized plas- 
mas. The differences between the conclusions of Gary & 
Feldman and Levinson & Eichler (1992) result from differ- 
ent magnitudes of the assumed magnetic contribution to to- 
tal pressures: Whereas Gary & Feldman assume comparable 
magnetic and gas pressures as is appropriate for the solar 
wind, Levinson & Eichler (1992) consider a weakly magne- 
tized plasma. 

Levinson & Eichler (1992) conclude that if the mag- 
netic pressure is small compared to the gas pressure, heat- 
flux inhibition along field lines is effective. In any case, both 
treatments neglect geometric details of the plasma processes 
involved. In particular, the assumption of a BGK scatter- 
ing operator used by Levinson & Eichler 1992 (followed by 
Pistinner et al. 1996) to calculate the particle distribution 
function, which neglects the pitch angle dependence of the 
scattering, was not explicitly justified. The purpose of this 
paper is to quantitatively assess the validity of the BGK 
model in view. We show that all pitch angles scatter at 
roughly comparable rates, so that the neglect of the pitch 
angles dependence (made by Gary & Feldman 1977; Levin- 
son & Eichler 1992, Pistinner et al. 1996) is quantitatively 
valid. 

The structure of this paper is: In § |^ we consider the gov- 
erning equations for collision dominated and whistler domi- 
nated plasma. In § ^ we obtain the formal expression for the 
electron distribution function first in the absence of whistlers 
and then in the presence of quasilinear situated whistlers. In 
§ H we obtain the critical condition for whistlers excitation 
by a collision dominated plasma. We then discuss the consis- 
tency of the small scale plasma turbulence theory, and show 
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that off-axis whistlers are required for significant heat flux 
inhibition and that they are excited by the required amount. 
The expression for the heat flux vector is considered in § ^| 
We summarize our conclusions in S 0. 



2 GOVERNING EQUATIONS 

This section outlines the governing equations that lead to 
the heat flux. The discussion starts with collision dominated 
plasma theory. It concentrates on the assumptions made in 
deriving the inhibited heat flux, in particular, the assump- 
tion of steady state, which raises the question of kinetic sta- 
bility of the plasma. 

Kinetic plasma instabilities can lead to a plasma per- 
vaded by stochastic electromagnetic fields. These fields may 
scatter particles faster than binary collisions. To account for 
a kinetic instability saturation, a set of governing equations 
that describe this phenomena is required. 

An electron-ion plasma, irrespective of whether the 
plasma is collision or collective phenomena dominated, can 
be described by two Boltzmann or Fokker-Planck equations 
(BFPE). For sake of brevity we shall assume that the solu- 
tion to the ion BFPE is given, and consider only the solution 
to the electron BFPE. The electron BFPE (e.g. Blandford 
k Eichler 1987) reads: 



^#7 + 



( 21 



+ 

d 1 



v{l~fi A ) dlnB Of 
2 din dfi 



(i) 

L -y~ v °j- 

dn 2 Sfj ' 

where /, fi, v, B, E® c , v, l\\ are the electrons distribution 
function, the cosine of the pitch angle, the electron speed, 
the magnetic field magnitude, the self consistent DC elec- 
tric field Q the electron scattering rate (collision frequency) 

and the affine length along the magnetic field lines respec- 
tively. 

Several assumptions have been made in writing eq. |l|: 

(i) steady state 

(ii) Lorentzian plasma, i.e. energy exchange between par- 
ticles is ignored 

(iii) curvature drifts are ignored 

(iv) phase velocity of the whister waves is ignored. 

These assumptions are justified in the context of clusters 
of galaxies, where the large scale hydrodynamical timescale 
and the collisional timescale are both long compared to the 
plasma timescales, the scale length of the magnetic field is 
enormous compared to particle gyroradii, and the whistler 
phase velocity is small compared to the electron thermal 
velocity. For completeness we elaborate briefly on each of 
these assumptions in § ^| 

Classical astrophysical plasmas [J] have traditionally 

* Note that Blandford & Eichler (1987) do not consider the self 

consistent DC electric field explicitly. In the case of cosmic rays 

the relative velocity of electrons compared to that of the ions is 

small in comparison to the velocity of the plasma. 

1 Note that the scattering rate of the electrons can be due to 

collisions with ion or collective phenomena (or both of them and 

it has not yet been specified). 

T With the exception of cosmic ray theory. 



been assumed to be thermal i.e. collision dominated. Thus, 
the Spitzer & Harm (1953) results are traditionally applied 
to model them. With eq. ^ one readily reproduces the Spitzer 
& Harm (1953) results (within a factor of five; cf. Table III 
and eq. 36 in Spitzer & Harm as energy exchange is ignored 
in eq. 1) along a single magnetic field line. Toward that goal 
one only need assume that Coulomb forces (collision domi- 
nated plasma) yield: 



V = Vcoll, 

where 

- _ v . __ 3 
and 

vtn = v C ou{vth) = hi(A coui ) 

V = 7T7 Vth = 



(2) 
(3) 

(4) 



with T, ks, rile, oJp, Nd and ln(A cou () the gas tempera- 
ture, the Boltzmann constant, the electron mass, the elec- 
tron plasma frequency, the Debye number, and the Coulomb 
logarithm respectively . 

To obtain significant electron heat-flux inhibition by 
whistlers, whistler-electron scattering must be as fast as 
ion-electron scattering. Thus, the time scale for establish- 
ing steady state on microscopic scales is at least as fast as 
that when whistlers are absent, and we may assume that a 
modified steady state is formed. Under such a situation eq. jj] 
describes the evolution of the electron distribution function 
but, with the following important modification, 



— 3 , - 



(•>) 

where ^14,= v w /uth, and v w is the electron whistler scat- 
tering rate. A few more implicit assumptions are made in 
writing eq. Q: The equation is written in a frame in which 
the whistlers' electric field vanishes. We assume that the 
magnitude of the DC field is not big enough to lead to in- 
elastic scattering events. The frame in which the scattering 
is elastic moves with the whistler phase velocity, which has 
been neglected in comparison to the velocity of the electrons 
exciting the whistlers. 

To obtain v w one has to choose an approximated non- 
linear plasma turbulence theory. In this paper we use quasi- 
linear theory. There are several formulation of quasilinear 
theory, and we chose to use the semi-classical formulation 
given by Melrose (1980). From his work, after some algebra, 
we obtain PJ 

v w = Q(n)u w>s + [1 — Q(ri)](v w , s + U w , a ) = 
Vw,s + [1 — Q(y)]{Vw,a)- 



(6) 



where 



6( M ) = 



n < 
n > 



(7) 



3 The details of the procedure that lead to the modifications to 
eq. |l| in which eq. ^ replaces eq. ^| are given by Melrose (1980). 
" Melrose 1980; VII, p. 20, f. 7.56; non-relativistic electrons 
assumed 
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1-^ 



W(k r 



Hlxl" ' 



(8) 



and ^(fe,^), f^e, fc, c, c, Xi Pi are the whistler spectrum, 
the electron gyro- frequency, the whistler wave number, the 
electron charge, the speed of light, the cosine of the angle be- 
tween the whistler propagation direction to the field line di- 
rection, and the electron momentum respectively. Although 
the above expressions can be written more compactly, (Mel- 
rose 1980), we write them in this way to reflect that v w%a , 
has contribution from asymmetric wave emission while u W)S 

has contribution from symmetric wave emission, [j] 

Equation ^ demonstrates that v is a function of 
W(k,x)- Thus, one requires an equation that accounts for 
the evolution of the energy density in the waves (or the 
waves spectrum) . Under the assumption of steady state the 
equation governing the evolution of the whistlers spectrum 



dW 



(9) 



where 7™ is the quasilinear growth rate (which is essen- 
tially the linear growth rate expanded under the assump- 
tion £ = 7™ /uj w <C 1), I™ is the non-linear damping rate 
and v g is the whistlers group velocity. In the framework 
of quasilinear theory the explicit expression for 7™ (Melrose 

1980) [fij reads: 



{p 3 v (yS£ 



• dp. 



+ 



V) 

i \ P dp 



)(>-**)}. 



where 



PR = 



m e Q, e 



(10) 



(11) 



is the resonant electron momentum and rm and va are the 
ion mass and the Alfven velocity respectively. 

We have a set of coupled functional equations that gov- 
ern the behavior of theplasma. These equations are en. hi for 
/, (with eq. | and eq. | ) and eq. | for W ( with eq. g). We 
emphasize that eq. ji] and eq. ^ will retain their form in any 
plasma turbulence theory, and that the turbulence theory 
determines eq. ^| eq. ^, and eq. |l^. With this set of coupled 
functional equations the distribution function can be deter- 
mined and the heat flux can be calculated self-consistently 

B 



II The goal of this paper is to establish that u WtS 2> 1 which 
is the implicit assumption made by Levinson &; Eichler. If only 
& w ,a 3> 1 then no heat flux inhibition is obtained. 

calculated for an arbitrary mode by second order perturbation 
theory (Melrose 1980); and explicitly for whistlers by Levinson 
1992 

tt Melrose 1980; VII p. 18 f. 7.54 

W The functional dependence is implicit: / is a function of u 
which is a functional of W which is a function of 7™ which is a 
functional of / etc. 



3 SOLUTION OF THE ELECTRONS BFPE 

We consider an approximate solution to eq. |l| under the fol- 
lowing assumptions: 

(i) We assume that a Knudsen expansion is valid, even in 
the absence of whistler waves. 

(ii) We assume that the net current vanishes. 

The first assumption implies that the distribution function 
can be expanded into a series with / = f° + f and that 
/° S> f 1 2> f 2 such that only the first order correction can 
be retained. 



3.1 The Knudsen Expansion 

The solution obtained by a Knudsen expansion can be vali- 
dated a posteriori, and this solution procedure is not sensi- 
tive to whistlers' absence or presence in the plasma. A lit- 
tle algebra shows that the solution to zeroth order is given 
by an isotropic distribution function which we take to be a 
Maxwell-Boltzmann distribution. With this assumption the 
first order correction (in dimensionless units) reads: 



th 



df _ 0. 

Oil V 

where 

(MB 



df 



eEi 



df 



dh 



m e vthVth.v 



dv 



f 

coll 



— A, 



: n e 7T 2 
th d\n{T) A 

Sill 



{m e vth) 



(12) 



(13) 



are the Maxwell-Boltzmann distribution, the Knudsen num- 
ber, and A* h = v t h/v t h is a thermal electron collision mean 
free path. Note that only f MB is still dimensional but all 
other quantities are dimensionless. To eliminate E DC from 
eq. ^| we need to relate it to the spatial variability of f MB . 
Toward that goal one uses the current free condition. 



3.2 The Current Free Distribution Function 

In astrophysical conditions, thermoelectric fields would 
cause the plasma to settle down to a zero current state. 
The current free condition reads: 



j\\ =2tvJ dvv 3 J dfifif = 
2n Jdvv 3 fd^ |£=0. 



(14) 



Using eq. ^| in eq. ^ and substituting the result into eq. ^ 
we find after some algebra that the current free distribution 
function is: 



21 

dfi 



v th 
— e 
V 



-2 



f 



(15) 



where 

Y 2 = f 
Jo 



dvv e e~* Z(v) 

yl = JT dvv A e-** Z(v) (16) 

and the ratio of Y 2 /Y 1 should be evaluated once the value 
of v has been determined. For a collision dominated plasma 
the value of v is given by eq. ^. Thus, the value of Y 2 /Y 1 
can be evaluated. However, if the plasma is whistler domi- 
nated, / is a functional of v, and not a simple function of it. 
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Substituting v = v 3 into eq. |l6|, one finds Y 2 /Y 1 —4 and 
using this in eq. one obtains 



f = f MB [1+^ 
f = .f MB [l + ^ th 



(17) 



This solution is valid provided that e th = e coU (v th ) < 2 x 
1CP 2 (due the velocity dependence of the collisions mean free 
path; Gray & Killkenny 1980). Note that if e th = 2 x 1(T 2 
then for v = 2 => e co " = 0.32 ^ 1, for v = 4 the underlying 
assumptions of the expansion do not hold. 



4 SOLUTION OF THE QUASILINEAR 
EQUATIONS 

The aim of this section is to a) calculate when eq. [l?] gives 
rise to whistlers, b) to find the angular dependence of the 
whistler spectrum, and c) with this whistler spectrum to 
estimate the modified steady state distribution function (to 
which the collision dominated plasma distribution function 
evolves), from which follows the heat flux inhibition factor. 
Consistency checks of the a priori assumptions made are 
carried out as we progress. To obtain these goals we consider 
approximate solutions to the equations given in § ^ under 
the assumption that whistlers pervade the plasma. However, 
before we do so we need to establish when this situation 
would occur and what the qualitative nature of it is. 



4.1 The Whistler Quasilinear Growth Rate 

The formal expression for the distribution function has been 
found via the Knudsen expansion and is given by eq. |l5[ 
To calculate when whistlers are present in the plasma we 
need to use eq. [if] ^| in eq. [lo| It proves mathematically 
convenient to formulate all expressions in terms of feii = k\ 
and x as the independent variables rather than k and \. 
After some algebra we obtain the following schematic form 
for the quasilinear growth rate: 

7 » = ^n e fc%P(k hX )[e th p e Q(k hX )p-\k hX )-l](lZ) 
where 

P(fc„ , X ) = 2^1 ft dtil - ^e~^ = 



Q{k hX ) = J d/i(l-M 2 )C 4 e" 5 

(i+x 2 )(i+Cnr)+(i-x) 



~4 ( ~1 
V I V 



1-H 



off 



+ 4 



(19) 



l.'r 



1 

fcy/j 

k = krf 

87rF e ~ la 



8ttP 
B 1 



SS One can derive the plasma dispersion relation and consider the 
growth rate for the whistler mode, under the assumption that the 
growth rate is much smaller then the oscillation frequency of the 
waves. This is tantamount to the use of the quasilinear growth 
rate. 



and P e is the electrons pressure, and P the total gas pres- 
sure. We have chosen the magnetic field pointing from hot 
to cold. For this choice whistler propagation is parallel to 
field lines, in the coldward direction. Thus, for wave growth 
we must have X > ( or > 0)- This point is taken into 
account in the expression written above. The terms Q{k\\, x) 
and P(k\\ , x) m ec l- represent a growth term due to the jM- 
term, and the damping term due to the ^ term in eq. |l0| 
respectively. 



4.2 Stability of the Collision Dominated Solution 

Here we show that under a wide range of conditions, a dis- 
tribution function, were it governed only by collisions, would 
be unstable to whistler waves which, of course, modify the 
distribution function. 

Using eq. 18 we first establish when a collision dom- 
inated plasma becomes unstable. Using Y 2 /Y 1 — 4 and 
u w ' s = u w ' a = in eq. |l|, we have: 



7£H = 2 £«.0r 1 (i + x a )fcu 



e th /3 e v 4 (v 2 - 4) 



/^(i-mV 



(20) 



where the subscript coll denotes whistler growth in a dis- 
tribution function stabilized by binary collisions only, not 
by the whistlers themselves. Some rough estimates may be 
obtained from eq. The first is that, provided e th fi e > 1, 
we have 7™ > 0, in a narrow cone along the axis. The width 
angle of this cone is roughly estimated to be /%r 1 ( e *' l )~ 1 < 1- 
To obtain the quantitative conditions we note that the inte- 
gral in eq. ^ can be written in terms of tabulated functions, 



(fc||) 4 e 



+ 



where 
erfc(x) = 1 



¥(*n 

(fe||) 7 er/cp|| 



— t 

dte 




(21) 



(22) 



Marginal stability of the plasma is found numerically in 
e e t h,ku,x space at 



(23) 



e th f3 e = 0.0001 ± 10~ 5 
jZuik = °- 22 ± 10_3 ' X = 1 ± 10- 4 ) = ± 10- 1 
The validity of the latter solution is subject to the condition 



0.00054 < PlxWe 1 < 1 
k 2 

0.00054 < < 1, 

where 

u w = n e Pp: 1 \ x \, 



(24) 



(25) 



which is the range in which the whistler mode exists. A 
numerical study of the properties of eq. ^ shows that there 
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is a distinct range of k\\ for which y™ a u > 0. This value is 
bounded between 0.22 < fcy < 0.7 for 0.0001 < e th e < 
1000. For fen = 0.22, one finds 



180 > P > 0.01. 

For fey = 0.7 instead of 0.22 we find that: 
1800 > p > 1. 



(26) 



(27) 



The underlying physical interpretation of this range of P is 
as follows: If P is too large, the waves that resonate with 
thermal electrons can cyclotron resonate with ions as well, 
which complicates things beyond the scope of this paper. 
For P too small, the phase velocity of the whistlers exceeds 
the thermal electron velocity by a sufficient margin that the 
resonant particles are far out on the tail of the Maxwellian 
distribution function, and the excitation is too weak to be 
interesting. In any case, condition eq. ^is met in clusters of 
galaxies, and condition [26] at the ISM. The angular depen- 
dence of 7 on eq. ^l] on x IS roughly linear with a maximum 
at X = 1; the slope depends on e th p e . Thus, at \ < 1j con- 
dition eq. |27] or condition eq. |26| are somewhat relaxed by a 
factor of order unity. 

The critical condition for whistlers present in a collision 
dominated hydrogen plasma (cf. the last line of eq. |l9| ) is 
found from eq. to be : 

e th p > 2 x 10~ 4 , (28) 

provided that P is in the range given by eqs. 

There are a few more conditions for the validity of 
eq. One of them is 

2 



€ th < 10" 2 



p > 10" 



(29) 



and it stems from the validity of the Knudsen expansion. 
The latter constraint is so weak that it can be ignored. 

The last condition to be considered is merely a formal 
one and it stems from the validity of eq. [lo] namely: 



^ s max(7^ oii ) ^ l 



(30) 



This condition is rather strong and it requires that e /3 e < 
10. However, even if this is not the case if when whistlers 
dominate the plasma, the value of v WjS is sufficiently high 
and £ can be maintained much smaller than unity. 

The conclusion of this section is that any hydrodynam- 
ical model in which P > 0.01 and 10~ 4 < e th < 10~ 2 gives 
rise to whistlers due to an excess of perpendicular momen- 
tum in the hotward velocity hemisphere. The unstable on- 
axis whistlers propagate in the cold ward direction. 

Note that the hotward part of the distribution function, 
which has excess perpendicular momentum is not the part 
that carries the heat. The coldward particles, which carry 
the heat from hot to cold, have an excess of parallel mo- 
mentum, and are stable to all on-axis whistlers 

Thus, the heat carrying particles cannot resonantly ex- 
cite on-axis whistlers that propagate parallel to the field 
lines and on-axis whistlers that propagate parallel to the 
field lines do not significantly inhibit heat flux. Whistlers 
that propagate off the axis, on the other hand, can inhibit 
heat flux. These off-axis waves have right-handed elliptical 
polarity which can be represented as a superposition of left- 
handed and right-handed circular polarity. The left-handed 



circular polarity allows for whistler emission in perpendic- 
ular momentum deficient regions in phase space. Although 
at marginal stability of a collision dominated plasma only 
on-axis whistlers are excited, it is incorrect to proceed un- 
der the assumption that these on-axis whistlers can lead to 
heat flux inhibition. Rather, we argue, any on-axis waves 
merely modify the distribution function in such a way that 
the marginal stability is transferred to off-axis waves. 



4.3 Condition for Whistler Dominated Plasma 
and Energy Density of Whistlers 

The questions are related as to when the plasma would 
become whistler dominated and how much energy density 
would be stored in the whistlers, and we consider them 
together. As was pointed out above, a whistler dominated 
plasma does not necessarily imply heat-flux inhibition. The 
plasma becomes whistler dominated when 



> 1, 



(31) 



whereas the whistler dominance leads to heat flux inhibition 
when 



> 1. 



(32) 



Using eq. |3| and the first line of eq. g in eq. [B2| we find 
that heat flux inhibition is obtained provided that: 



£ R 2 9 Nd 



(^)/ 1 1 _^(^.x)(i-x 2 )^»i. 



(33) 



Writing the spectrum in the following general way: 

W = E i K i (k ]] )S\ X ), (34) 

where K l are arbitrary functions of ku, and H 8 (x) are arbi- 
trary functions of \ we have: 



'3? 



(35) 



Thus, if £ is sufficiently big we have heat flux inhibi- 
tion. Assuming that this is the case, we write W(kR,x) = 
Const W(k) and assume Const is of order unity. We can 
now invert the above ratio and obtain the constraint on the 
energy density in the waves which reads: 



4.nU'(k n )k n 




(36) 



W(k) = N-^lpiLogiA^o) . 

The energy density in the whistlers is a function of the in- 
hibition factor, and inversely proportional to it. Since the 
inverse of the Debye number is tiny, 10 -16 for the ISM and 
10 -23 for the ICM, inhibition factors can be extremely large. 
Thus, even an inhibition of the heat flux by a factor of 10 10 
leaves the energy density in the whistlers small compared to 
the energy density in the magnetic field. Therefore, one can 
quite safely assume that the energy density stored in the 
whistlers is invariably much less than the energy density of 
the background magnetic field. 
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4.4 Whistler Dominated Plasma 

We have shown that once the condition is satisfied the 
plasma becomes whistler dominated. We now proceed under 
the assumption that eq. ^ is satisfied by a large margin, 
i.e. that 



e th p > 2 x 1(T 4 . 



(37) 



Under this conditions whistlers grow as 7™;; > 0, for at least 
the same mode, and we need to solve the wave equation eq. ^| 
Thus, we need to determine W, and v self-consistently. We 
consider quasilinear saturation i.e 7™ = 0, for the most un- 
stable mode. In particular we need to know whether quasi- 
linear saturations sets in before or after eq. |32| becomes rela- 
vant, as only the latter leads to heat flux inhibition. Eq. |3lJ 
though implying a whistler dominated plasma, merely im- 
plies a reduction of the heat-flux only by a factor of two. 

Advection by waves is negligible (cf. § |^) and the wave 
equation eq. ^ with eq. [73] reduces to 

j w =T w . (38) 

Quasilinear saturation is defined as the vanishing of 7™ and 
hence of T w . If the wave spectrum achieves quasilinear sat- 
uration to a good approximation, i.e. for every k\\ there is a 
critical angle Xcnt such that: 



7™[fc|hXcWt(fc||)] = 0, 



(39) 



and for all other angles the growth rate 7™ is negative. Thus, 



W OC 8(X~ Xcrit). 



(40) 



The picture is that the wave spectrum is highly peaked near 
Xcnt, where the condition of marginal stability is met. Once 
this has been established, the condition £ <C 1 (cf. eq. |30[ ) is 
a priori fulfilled. 

The assumption of quasilinear saturation is exact if 
Xcrit is unique and does not depend on ku . Otherwise the 
assumption F w = is formally invalid (Mode coupling does 
not exist if all the waves propagate at the same angle) . How- 
ever, the non-linear Landau damping rate F LD is of order 
P~ 2 compared to the mode coupling term, and it can be 
shown (Levinson 1992; cf. |l^) that if the parent wave and 
the resulting waves have energies of the same order, then 
the non-linear damping is less than u w ' s /v coll r] (r/ is defined 
in eq. [70] ). Thus, it is of order the suppression factor times 
r\ and can quite safely be ignored for the suppression factors 
that we are interested in. 

Using eq. [iS] in eq. ^ we obtain: 

e^PeQik'XcriJP^iku, Xcrit) -1 = (41) 

We now seek solutions of eq. ^ which will yield a value of 
Xcnt by exploiting the fact that Q(k\\ , Xcnt)P~ 1 (k\\ , Xcnt) 
has a single maximum of x a t Xcnt- To start this proce- 
dure we evaluate Q, which is defined in eq. In particu- 
lar we note the Q dependence on the ratios v w ' a jv col \ and 

. „~,-3 



we find: 



at it 



Using eq. 34 in eq. H, and u c 



(42) 



= H 



(fc ||M )3 



£iif(fc||)<Si 



where 



3iV D /3 e 



41n(A coui ) V>"e \ c I ^ 

5, = J o l d X ^(l-x) 2 . 



(43) 



We consider now eq. ^ which gives the relevant case. The 
elimination of the case given by eq. ^l| is considered further 
below. 

A whistler dominated plasma under the condition eq. ^ 
implies that #AiK(k\\) > 1 and HS,if(fc||) > 1 for the 
dominant i. Substituting eq. Wa into eqJlSl we find 



(1 + X 2 ) + (1- X ) 



( fc ||M) fc|| 
S i K i (k,,)(A i +S i ) 



M) 2 

2 EiiCi (£11 ).Aj 



(44) 



Using the last equation in eq. ^ we obtain: 



th o 



G"k 



f S,K( fc|| )A Xc „ t( i- Xe „ t) - 



(45) 



where 



G n (fc,|) = J-e^ 

K11 ./() 



d^l-fl 2 ) 



(fc||M) 2 



y2 

Y 1 



,(46) 



is a rational function of fcii , and 



r s _ 1 l-s 
G ~2 II 



(47) 



Substituting eq. |^| into eq. ^ one finds 

i^Xcri 



2KE i if(fc|,)(S i +^ i ) 



1 I ~ i, 



i+x: 



If we proceed under the assumption that only a single value 
of i contributes to the sum, we find: 



th o 



1 



2HK(k 



)S(1+K) i^cri 



+ 2- 



Y' 



k'i ' ~ Y 1 



(49) 



(50) 



where 

We can now obtain the critical angle for whistlers propa- 
gation at marginal quasilinear stability. Note that the func- 
tion given in the parenthesis might have a local maximum 
in Xcrit- One finds this maximum from the condition: 



1 , (1-XeWtr I 

^ + 1+x 2 + 

it(l-4XcWt+3x 2 wt ) 2xf„ t (l-Xcrit) 2 



1+X 







(51) 



Thus, the location of the maximum for the value of Xcrit de- 
pends on the value 1Z only. Using eq. |4o| in eq. ^ to estimate 
the value of 1Z we find: 



1Z 1 



^Xcrit 
(I -Xcrit) 2 ' 



(52) 
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The cubic equation eq. ^l] has two complex roots and one 
real root. Substituting into the expression for the real root 
the value of 1Z from eq. we find an equation for x^nt- 
We solve this equation numerically, and find that 



n = 4.8, 

and 

Xcru = 0.414214 



« 64°. 



(53) 



(54) 



We now substitute the values of Xcrit and 1Z into eq. |49| and 
find: 



th rt 



Y 2 



Y 1 I 2HK{k«)S 



(55) 



One can now solve this equation for the value of Ki^(fe|i)>S. 
Bearing in mind that only one i contributes to the sum in 
eq. ^ and substituting the value of ~^.K(k^)S into the last 
line of eq. H2| we find: 



3e th f3 e 
fj,v 



(Vflf 



Y 2 



(56) 



and with the numerical value of 1Z we have from the first 
line of eq. 



Ue th p e 



(« M ) 2 + 2 - 



(57) 



which implies that the heat flux along field lines is roughly 
proportional to fl^ 1 , and can be determined from observable 
parameters. The solution procedure continues by finding the 
value of the ratio Y 2 /Y 1 . This problem can be considered 
separately, and we do so in § Thus, with eqs. ^ and 
eq. [l5| we have: 



K 1 (i 



8£ 



3+U[l-B(fi)] 

(m«) 2 +2-| 



(58) 



and from § [u] we find Y 2 /Y 1 
function reads: 



0.18. Thus, the distribution 



3+14[1-G(m)] 



(59) 



AT 1 (v 2 - 0.18) [(/i«) 2 + 1.82] 

and a self consistent distribution function at quasilinear sat- 
uration of off-axis waves has been found. 

Although in a collision dominated plasma whistlers are 
least stable along the axis we now show that at quasilinear 
(QL) saturation this is no longer the case: If the waves at QL 
saturation would have propagated in a narrow cone along the 
axis this would have implied that HAiKi(k\\) ^ 1 f° r some i 
and that Si <C Ai for every i. Substituting this assumption 
into eq. ^ we find: 



th q 



Xcrit 



G 5 (fc||)fcf 
KEiif(fe||)A 



+ G 8 (fe,|)^ 



rit(l — Xcrit) 2 
l + y 2 . 



1(60) 



Now G 5 and G 8 are functions of the same order of magni- 
tude. The function in the parenthesis has a local maximum 
at Xcrit- Provided we neglect the term of order 1/H the loca- 
tion of the maximum does not depend on k\\ or e th /3 and it is 
found at Xcrit ~ 0.3, in contradiction with our assumption 
that x^u ~ 1 that Si is of order Ai. Since the corrections to 



the location of the maximum would only be of order 1/N, we 
conclude that it is very difficult to support a whistler domi- 
nated plasma that does not generate off-axis waves without 
imposing strict fine tuning. 



5 THE HEAT-FLUX ALONG FIELD LINES 

The heat flux vector is, after integration by parts, 



911 



ilrr ' I d/J,(l — [l 2 )^-. 



When eq. ^9] is used in eq. [il] one gets: 

5,1 = - 7 r-im e u e 3 n e e th (/3 e e*' 1 )- 1 
/* dvv 7 e-° 2 (v 2 - 0.18) £ d/i(l - ^ 2 )fi 
3+i4 [1 -e(M)] ^" 1 M a + 1.82] _1 . 
Carrying the integration over the angles yields 

1\\ = f ^-mevln^iP^y 1 



(61) 



(62) 



f °° dvv 3 e-* 2 



(63) 



Carrying the integral over the velocities we find: 
qy = 0.13ir- : 5m e v*n e e th (l3 e e th )- 1 

T(3.5) - 0.18r(2.5) + 0.083(-2(1.82) + (1.82) 2 + 1.82 3 



6e 182 Gf'° I 1.82 




1,1,4,0 
n e e th (p e e 

_L q _i 

-0.657T 2m e vin e p e , 



-0.65ir-2m eV 3 ne e th (f3ee th )- 1 



where 



<- T n 



ai, . 

6i,. 



2^7 9l asz n* , =ji+1 r(a i+s )nf ii+1 r(i-b i - s ) 

is the Meijer G function ( Gradshteyn & Ryzhik 1980; p. 897, 
f. 4), and (f> denotes an empty set. The numerical evaluation 
of G was carried out by the use of Mathematica. 

The last equation shows the heat flux strongly sup- 
pressed compared to the Spitzer value when the scattering 
rate by off-axis waves exceeds the electron-ion collision fre- 
quency (Levinson & Eichler 1992). Previous literature for 
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heat flux inhibiting astro-physical environments used a term 
/ s " pp , the so called suppression factor, that multiplied the 
collision dominated plasma heat flux i.e. 



9|| = -f supp 6^€ tn n-^m e vtn e , (66) 

where we used eq. [l?] in eq. |6l] to obtain eq. ^ and the factor 
jsupp wag j n ^ roc j ucec i D y hand. The value of f supp obtained 
from eq. 64 and p6 is 



f supp = (64/65) x KTV^r 1 ~ 1 x 10" 



2 th n 



(67) 



If the distribution function were not to relax to quasilin- 
ear marginal stability then there would be less suppression. 
However, the non-linear damping terms have been shown to 
be much smaller than the quasilinear damping term (which 
in marginal stability is balanced by a growth term) and thus 
it can safely be argued that the suppression factor computed 
above is of general applicability (provided that f3 satisfies the 
required constraints). Since f supp has been calculated under 
the assumption eq. we suggest the following extrapola- 
tion: 



r 



1 + 50(e th /3) ■ 



(68) 



Note that the suppression factor has been calculated relative 
to a Lorentzian plasma. Taking into account that the Spizer- 
Harm heat flux is about a factor of five higher relative to a 
Lorentzian plasma, the suppression factor reads: 



r 



i 



1 + 250(e th /?) ' 



(69) 



6 DISCUSSION AND CONCLUSIONS 

Heat-flux inhibition along field lines obtains when off-axis 
waves are present in the plasma. We have found within 
quasilinear theory that whistler instability indeed gener- 
ates off-axis whistlers. Thus, heat flux inhibition along field 
lines by whistlers should be a common phenomenon in a 
weakly magnetized plasma. The basic argument for off-axis 
whistlers is, in qualitative terms, that if the waves were ex- 
cited only in a narrow cone along the axis, then the distri- 
bution function would evolve to a state in which the growth 
rate of off-axis waves would generally exceed that of the on- 
axis ones. 

One of our assumptions is that the energy density in 
the magnetic field of the whistler will not exceed the energy 
density of the background magnetic field. We find that this is 
generally the case by a large margin in tenuous astrophysical 
plasmas. 

While Gary & Feldman (1977) find whistlers inefficient 
for inhibition of heat flux in the solar wind, f aupp ^ \ } Levin- 
son & Eichler (1992) find a suppression factor of 1CT 2 < 
fsupp < 10 -i ) in the interste ii ar medium (ISM), sufficient 

to explain ISM observations. Pistinner et al. (1996) find 
1 > f supp > 10" 4 in the intracluster medium, ICM, which is 
enough to resolve many of the observational problems there. 
In these papers, the level of wave turbulence was calculated 
by assuming a balance between linear growth and non-linear 
wave removal, and it was noted that the results are suspect 
when the heat flux so found was below that needed to drive 
the whistlers beyond quasilinear saturation. Here, on the 



other hand, we have used the criterion of marginal stabil- 
ity of the least stable whistler wave to derive a relatively 
simple expression for the heat flux [eq. or, equivalently, 
the suppression factor by which classical heat conduction is 
suppressed of eq. pS[ . The heat flux that we derive here is 
at least that found in our previous papers inasmuch as of 
the geometry of wave - particle interactions has been con- 
sidered more carefully, i.e. only some particles can interact 
with some waves, as opposed to all particles with all waves 
via an oversimplified BGK operator. On the other hand, it 
will be argued in a subsequent paper that the asymmetry of 
the scattering in the forward and backward hemispheres is 
sufficient to allow for a non-resonant firehose instability of 
the magnetic field. This causes field line tangling, and fur- 
ther inhibits the heat flux. A careful evaluation of the heat 
flux as a function of plasma parameters will therefore be 
attempted only in the subsequent paper. 
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8 APPENDIX A 

Here we briefly outline the nature of the assumptions made 
in writing the BFPE eq. The assumption of steady state 
does not imply that the plasma does not evolve with time. 
Rather it implies that any microscopic process occurs on 
a time scale which is much faster than the time scale 
for macroscopic changes and thus, temporal macroscopic 
changes may be ignored to first order. 

The second assumption implies that energy exchange 
in a collision is negligibly small and may therefore be ig- 
nored. If the electrons collide only with ions it implies that 
the ions have infinite mass. Thus, electron-electron collisions 
have been ignored. If the plasma is collision dominated this 
assumption modifies the collision dominated plasma distri- 
bution function by up to a factor of five. A plasma dom- 
inated by collective phenomena which leads only to pitch 
angle scattering (the case considered in this paper) is tan- 
tamount to the above approximation. 

The third assumption is related to the ratio of the gyro- 
radius and the large scale magnetic field variation, the so 
called correlation length (not to be confused with the co- 
herence length). It implies that this ratio can be effectively 
set to zero. There is a subtle issue involved. If the plasma 
is pervaded by collective magnetic phenomena (electromag- 
netic waves) that have tl ~ Ib it seems that the approxima- 
tion which led to eq. |l| is not valid. However, this presents 
a technical problem that can be circumvented by a stan- 
dard procedure which allows one to include the small scale 
magnetic fluctuation into v (cf. Melrose 1980) . Observations 
suggest that the large scale magnetic field variation i.e. the 
magnetic field correlation length, is comparable to the elec- 
tron mean free path. Thus, we may define: 
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n 



L B 



n 



where A e is the electron mean free path, and Lb the mag- 
netic field correlation length. Typical values of rj are: 

(i) T) = 10" 12 ISM 

(ii) r) = lCT 16 ICM. 

Assuming that fj w ? 
approximation" . 



(70) Consider now the ratio of the quasilinear damping term to 
the non-linear damping term. 



one gets the "adiabatic drift 



9 APPENDIX B 

Consider the advection of whistlers. The magnitude of the 
quasilinear damping term Q, e f3~ 1 kk\\P{k, x) ( ec l- 0) com- 
pared to the advection term in eq. H is: 



(71) 



The entire whistler emission process is driven by gradients 
of the distribution function. Thus, the shortest length scale 
for the variation of the spectrum is determined by the length 
scale for the variability of the distribution function and we 
have: 



XV g dln(W) 
2 Sin 



Q, 



M _1 xlxl ( 
n.k0t 1 x \x\v(? k 



dln(W) 

Sill 



On the axis: 



10" 



(72) 



(73) 



10 APPENDIX C 

Levinson (1992) shows that non-linear Landau damping is 
given by: 



/327rm e u 



■ J dxidkiuji 



vy(fci.xi)^™ „ 2xi840(fc N 2 v a )* 



1 : |Sr? « 10- 8 « 1 
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(79) 



Y 1 



obtained. Although in principle we deal with an integral 
equation, it can be reduced by the procedure outlined below. 
Let 



a = 2 



Y 1 ' 



then substituting, eqs Hand [57] into eq. |l6| we find, 
Z(v) = (2e th p e )~ 1 £i- 3 



(80) 



(81) 



One can now proceed and calculate formally the value of Y 2 
and Y 1 . Substituting the expression for Z(v) into eq. [l6| one 
finds after some algebra 



f e th (5 e Y 2 



1 - 0.25[a + a 2 + 2e a G^° 



e G, 



0,1,1 
0,0,3,0 



(82) 



f e^PeY 1 = 0.5[1 - e a G 4 4 '° a 



0,1,1 
0,0,2,0 



Subtracting from 2 the formal ratio of Y 2 /Y 1 one obtains a 
transcendental equation for a. This equation reads: 



43TTl e f^|fc|| 2 I 
, , ,1V 1 Ml 

fc|| = fe]],! + fc|| )2 . 

Casting the last equation into a non-dimensional form we 
obtain: 



(74) 



a- 2 + 2 



1 - 0.25(a + a 2 + 2e a G 4 ]o a 



1843/3 C 



e &,' n a 



1 - e a Gt'° Q a 




0,1,1 
0,0,3,0 



0. 



(83) 



0.14&- 2 / dxidk 



(75) 



0,0,2,0 

Solving this equation numericaly, we find a = 1.82. 



The mode coupling term (Levinson 1992), after it is cast 
into non dimensional form reads: 



~,MC 



16 

Thus, 



— (^e^ 1 ) / do»rdfei(i-xi) ; 



(^) 5 



W(fci, X )(76) 



(77) 



Levinson (1992) evaluates the mode coupling term under the 
assumption that the parent waves and the daughter waves 
have energy of the same order and finds: 



r-«i J ' <(n^- 1 )*f^. 



(78) 
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